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Abstract 

The Yang Mills equations provide a classical mean field description of gauge 
fields. In view of developing a coherent description of the formation of the quark 
gluon plasma in high energetic nucleus-nucleus collisions we study pure gauge 
held dynamics in 3+1 dimensions. In collisions of wave packets, numerically 
simulated on a SU(2) gauge lattice, we study transverse and longitudinal en- 
ergy currents. For wave packets with different polarizations in color space, we 
observe a time delayed fragmentation after the collision resulting in a rapid ex- 
pansion into transverse directions. We call this phenomenon the "glue burst". 
An analysis of the Yang Mills equations reveals the explanation for this behav- 
ior. We point out that this effect could play a role in ultra-relativistic heavy-ion 
collisions. 



I. INTRODUCTION 

It is one of the most challenging topics in the theory of ultra-relativistic heavy-ion collisions to 
develop a coherent description of the formation of the quark gluon plasma. Many descriptions 
of the time evolution of such collisions have been developed in the past years [1—5]. Some of 
the most recent are the ultra-relativistic quantum molecular dynamics model (UrQMD) [1] and 
the parton cascade model [6]. While the UrQMD has been developed for collisions at center 
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of mass energies below 10 GeV/u, the parton cascade model is valid for energies at and above 
lOOGeV/u where the transition into the plasma phase is expected to occur. Both are based 
on the idea of a perturbative scattering of particles within transport models and describe in 
principle the evolution of a collision from the first contact of the nuclei throughout the high 
density phase until to the point where the last Hadrons are freezing out. These models however 
contain still a variety of problems. One of the problems concerns the description of the initial 
state of the colliding nuclei. The transport equations start from probability distributions of 
partons in the phase space. In reality however, the states of the nuclei are described through 
coherent parton wave functions. The incoherent parton description especially breaks down at 
exchanges of small transverse momenta. A few years ago, McLerran and Venugopalan proposed 
[7] that the proper solution of these difficulties is the perturbative expansion not around the 
empty QCD vacuum but around a vacuum of the mean color fields which accompany the 
quarks in the colliding nuclei. This idea motivates the development of a combination of the 
parton cascade model [6] with a coherent description of the initial states. This difficult goal 
can be accomplished in small steps through systematic studies. Recently, we have proposed 
a combination of a gauge lattice description for the soft color fields with a transport model 
for color charged particles [8]. Leaving out the collision terms first, this model then has been 
applied to simulate the collision of clouds of color charged particles accompanied by soft color 
fields in 3+1 dimensions [9]. The field energy distributions obtained for times shortly after the 
collision have shown transverse energy flows resulting from glue field scattering in the center 
of collision. For times larger than lfm/c glue field radiation seemed to be dominant. The 
sudden appearance of transverse energy flows during the overlap time of the nuclei was the 
motivation to study here the pure glue mean field dynamics leaving out the particles. The 
time evolution of colliding Yang Mills field wave packets has been studied a few years ago in 
1+1 dimensions [10]. These calculations have shown that wave packets of equal polarization in 
color space (Abelian case) do not interact whereas wave packets polarized into different color 
directions (non- Abelian case) decay after the collision into low frequency modes of high particle 
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multiplicity in the sense of Ref. [11,12]. This mechanism of a coupling between high frequency 
short wavelength modes and low frequency long wavelength modes in the Yang Mills equations 
has also been observed in Ref. [13]. 

In the present paper, we focus on the transverse dynamics and the coupling between longitudinal 
and transverse energy flows in collisions of Yang Mills fields. A study of the transverse dynamics 
requires simulations in at least 2+1 dimensions. Subsequently, we describe the method used to 
solve the Yang Mills equations and present results obtained from collisions simulated in 3+1 
dimensions on a SU(2) gauge lattice. These studies also reveal an interesting behavior of the 
time-evolution of non-Abelian gauge fields. 

II. TIME EVOLUTION ON THE GAUGE LATTICE 

In the Lie algebra LSU(2), we define the adjoint gauge fields *4. M (x) := A^(x)T c and the adjoint 
field strength tensor T^ v := F^ U T C . Einsteins sum convention has to be applied in the Euclidean 
metric for upper and lower color indices and in the Mincowski metric for upper and lower greek 
indices. The symbols T c with color index c = 1,2,3 denote the generators of LSU(2) obeying 
the commutation relations [T a ,T b }_ = if a b c T c and hence A tl {x),J- tw G LSU(2) for all xeR 4 . 
Here, we chose the representation T c = 1/2 a c with the Pauli matrices a c . Further below, we use 
also T° := |l 2 which is linearly independent from the generators T c . With these conventions, 
we denote the Yang Mills equations in the short form 

[2^1- = 0, (1) 

where is defined 

V, := d» - igA». (2) 

With this definition of the covariant derivative on the SU(2) mainfold, and with the def- 
initions S^x) := E£(x)T c , B»(x) := B^(x)T c (£^(x),B^(x) G LSU(2)) of the adjoint color 
electric and color magnetic field quantities, the Yang Mills equations Eq. (1) can be expressed 
in a form which is similar to the U(l) Maxwell equations in the vacuum 
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[D,£]_ = 0, (3) 

[D,B}_ = 0, (4) 

[Vx,£]_ = [V ,B}_, (5) 

[Vx,B]_ = [V ,S]. (6) 

With the condition [A(t,x),A(t,x')]_ = for all x, x' G R 3 for one arbitrary real time t, 
the equations (3) to (6) are indeed identical with the Maxwell equations. In this so called 
Abelian (or U(l)) case, one expects a linear behavior of the solution and it therefore provides 
an important test through comparison with the solution in the general non-Abelian case. The 
Yang Mills equations can be solved in an efficient manner on a gauge lattice in a Hamiltonian 
framework where we choose the temporal gauge A = 0. 

A lattice version of the continuum Yang Mills equations is constructed by expressing the color 
field amplitudes as elements of the corresponding Lie algebra, i.e. £g,k,Bx,k £ LSU(2) at each 
lattice site x. We define the following variables. 

Z4,i = exp{-iga t Ai(x)) = Z4 + , _, (7) 
Ux,u — Mx,kMx+k,iMx+k+i,-kMx+i,-i (8) 

In adjoint representation the color electric and color magnetic fields are expressed in terms of 
the above defined link variables Us,i and plaquette variables Us,m in the following way. 

ElJ = -Lu, Mj (9) 

®s,j = ej-fci {Us M - ui kl ). (10) 

The lattice constant in the spatial direction I is denoted by a\. As one can see from (6), the 
gauge field As,i is expressed in terms of the link variables Ugj e SU(2), which represent the 
parallel transport of a field amplitude from a site xeX to a neighboring site (x + l)eX in the 
direction /. We choose Us,% and £s,i as the basic dynamic field variables and numerically solve 
the following equations of motion. 
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Us,k{t) =iga>k £s,k(t) Z4,fc(i) 



£x,k(t) 



E{^,«(*)-z4,h(*) 



2ga 1 a 2 a 3 l=1 
Ul_ hl {t) (%_,,«(*) - z4_,,«(f))Z4-i,i(*) }■ 



(11) 



(12) 



III. NUMERICAL RESULTS 

First, we study the collision of plane wave packets which have initially constant amplitude in 
the transverse planes of a 3-torus lattice implying periodic boundary conditions. The size is 
chosen 20x20x2000 lattice points. We denote the lattice spacing in the space direction I by 
a\. For a collision, we initially arrange wave packets of a Gaussian shape in the longitudinal 
direction. The average momenta of the wave packets are (0, 0, ±/c 3 ) and their momentum spread 
is denoted AA; 3 . As a consequence of the lattice discretization these quantities are restricted to 
values between the maximal and minimal Fourier-momenta on the lattice: 

7, min _ 1 imax _ 1 /iq\ 

k™ in < A/c; < jfe, < fc[" aa: (14) 

The polarization in color space is defined by the unit vectors n£ for the left (L) and for the 

and right (R) moving wave packet respectively. 

We express the initial conditions through the gauge fields 

A c L = b lx n c L ^t,x ?> -Z L ), (15) 
A c R = 5nn c R <f)(-t,x 3 -Z R ) (16) 

where the scalar function <p(t, x 3 ) defines the initial wave packet 

(j)(t,x 3 ) =0o exp (-^Ak 2 3 (t + x 3 ) 2 ) cos (k~ 3 (t + x 3 )), (17) 



(2AA; 3 )/(v^<jA;3)(1 + e -^ Ahi > ) (18) 

The amplitude factor O is determined for a normalized wave packet when a = 100aia2. The 
parameter a is used to control the amplitude and allows to adjust the energy contained in a wave 



packet. In a particle interpretation it describes the total cross section per particle contained in 
the wave packet [10]. The second term in O is usually negligibly small for wave packets since 
k 3 /Ak 3 ^> 1. We chose the polarizations in color space ur = Hl = (0, 0, 1) for the Abelian case 
and n R = (0,0, 1), ul = (0, 1,0) for the non-Abelian case. Once the initial fields are mapped 
on the lattice, the time evolution of the collision starts from the superposed initial conditions 

41 = 41 R + 41 l, (20) 

respectively. A linear superposition of solutions obeys the Yang Mills equations only if these 
solutions have no overlap . Therefore, the initial separation AZ should be much larger than 
1/Ak 3 . 

The calculation contains four parameters. The relative color polarization which is parameterized 
by the angle 9q defined through n L ■ n R = cos(#c)- The average momentum of the wave 
packets k and their width Ak 3 and the coupling constant g which can be rewritten in terms 
of the parameter a as g' — g/y/cr by simultaneously rescaling the field A as A' = y/aA. 
Consequently, the system shows the same dynamics for different values of g and a, as long as 
the ratio g' = g/y/a is kept fixed. Multiplying the Eq. (32) in section 4 by a factor g, shows 
that each amplitude A^ absorbs a g such that no g is left over in Eq. (32). For /c 3 3> AA; 3 , 
the energy contained in one wave packet is essentially k 3 /a, i.e. for given k 3 the energy is 
determined through the amplitude and vice versa the amplitude is thus fixed by the energy. 
Once the amplitude is fixed, it makes sense to use the coupling constant as an additional 
parameter which determines the reaction dynamics of colliding wave packets. This is preferable 
also in view of studies in the future where normalized Dirac fields or color charged classical 
particles shall be included. In such cases g appears also in front of the source current of the 
inhomogeneous Yang Mills equations. We therefore keep the explicit denotation of g in the 
subsequent sections. 

In the following we present results from a simulated collision in the non-Abelian case with 
the parameters k 3 = n/(2a 3 ), Ak 3 = 7r/(100a), g — 1, a — 100a 2 = l.Ofm 2 . The simulation 
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was performed on a uniform lattice with constants a\ = 0,2 = as = a using a time step size 
At = a/10. If we would consider a collision of nuclei at low energies, a = 0.1 fm would 
be an appropriate choice for the above parameter settings. In particluar, A/c 3 = 7r/(100a) 
corresponds thus to a FWHM of 10 fm which is approximately the diameter of a 208 P6-nucleus. 
It is important to note that for the above chosen values, the calculation runs in the regime of 
weak coupling. Subsequently, we refer to the direction of the collision axis as the "longitudinal 
direction" and to directions perpendicular to the collision axis as the "transverse directions". 
Accordingly, we define the transverse and longitudinal energy densities of the color electric field 



Fig. 1 displays Wj> plotted over the ^-coordinate at various time steps t n as indicated on top 
of the curves. 
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FIG. 1. Transverse energy densities w^\t, £3) are plotted over the collision axis (longitudinal 



coordinate) for selected time steps t n . 
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At initial time to the distributions (to, xs) of the wave packets are completely filled resulting 
from the strong oscillations in the longitudinal direction according to = n/(2a 3 ). After 1500 
time steps the wave packets are colliding and have reached maximum overlap. The smooth 
white zone at the bottom results from a phase shift between the two superposed waves, i.e. 
the maxima of the wave packet (1) to not coincide with the maxima of the wave packet (2) at 
time step £1500- After the wave packets have passed through each other (at about £2500)1 a small 
white zone remains in the distribution w!^\x 3 ) of each and keeps continuously growing. At 
time step t 60 oo (t = 600a) the large fraction of the initial high frequency oscillations is reduced 
to a small remaining contribution visible on the surface of the distibutions in Fig. 1. The hight 
of the two receding humps decreases accordingly while the energy carried by each wave packet 
is constant in time. Almost all the energy which has originally been carried by short wavelength 
modes around k 3 has been transmitted into long wavelength modes which have filled up the 
vallies in the oscillating distribution w^\x 3 ). This behavior agrees qualitatively with results 
obtained in Ref. [10] where collisions of wave packets have been studied on a one dimensional 
gauge lattice and for times not larger t = 600a. At time step tsooo however, we observe that 
energy is partly transmitted back into high frequent modes before the wave packets start to 
decay around time step £10000 • At time step t^ooo, the energy distribution has expanded into 
longitudinal direction. The appearance of a circular polarization of the receding wave packtes 
could in princple lead to the same behavior before they decay. This possibility, however, is 
excluded by our numerical results which show - in the case of plane waves - no excitation of 
modes with x 2 -polarization for all times. This can be veryfied down to floating point precision 
but it has also a fully analytic explanation through Eq. (32) in section 4 for v = 2. 
This behavior does not appear in a collision of wave packets which are equally polarized in color 
space. In this case the shape of the two humps in the distribution Wj?\xs) at the time step 
£14000 is almost identical with the shape at time to- Deviations result from lattice dispersion. 
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FIG. 2. Longitudinal energy densities w L (t,x^) are plotted over the collision axis (longitudinal 
coordinate) for selected time steps t n 

In Fig. 2, the corresponding longitudinal energy densities w^\t,x 3 ) are displayed on a log- 
arithmic scale for selected time steps t n . We remember that the wave packets were initially 
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polarized into the transverse Xi-direction and consequently w L has to be zero as long as they 
propagate free. However, when the colliding wave packets of different color start overlapping, 
we observe an increasing longitudinal energy density in the overlap region around the center 
of collision at rr 3 = 100 fm. Fig. 2 clearly shows that Wj^\t,z) grows rather fast from time 
step tiooo until time step t^oo where it reaches a maximum and has grown by more than two 
orders of magnitude. For larger times the hump decreases again and practically disappears 
at £3000- Around teooo, however, the longitudinal energy density grows again at the positions 
of the receding wave packets. After 10000 time steps w^\z) has increased by five orders of 
magnitude. 

A comparison of Fig. 1 with Fig. 2 indicates that the total color electric field energy on 
transverse links of the lattice decreases at large times while the total color electric field energy 
on longitudinal links increases. In order to understand this behavior in detail, we explore the 
time dependence of these quantities for different values of the coupling constant. To be more 
precise, we define the transverse and longitudinal energy of the color electric fields as 



The integration in the definitions (23) and (24) is carried out over the whole lattice. 
Fig. 3 displays Wj, (t) over a time interval of 600a which corresponds to 6000 time steps. The 
energies are compared for different values of the coupling constant g. As already mentioned 
above, values of the coupling constant which differ from g — 1 can be scaled out with the 
amplitude of the wave functions. The amplitude of wave packets, however, determines the 
energy carried by a wave. Once the energy is fixed by describing a real physical system, it 
makes sense to use different values of the coupling constant. Here, we collide wave packets 
which are normalized to NiN 2 a 2 /a particles and refer to different values of g. 
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Fig. 3 
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FIG. 3. The transverse color electric field energy Wj, '(t) as a function of time is displayed for dif- 
ferent values of the coupling constant g. The time step width of each time step is At = 03/10 = 0.01 fm. 

Fig. 3 shows that for small coupling, Wj~ (t) remains unchanged through a long period of 
time after the collision which occurs in an interval around the time step £1500- F° r 9 — 4 the 
transverse color electric energy starts to decrease around time step £5000- A comparison with 
curves obtained for increasing values of g shows that the Wj< (t) begins to decrease at earlier 
times. For the largest values g = 10, 12, 14, the decrease begins in the overlap region of the 
wave packets. In these cases, strong oscillations occur in the overlap region. The magnetic 

( R) 

field energies Wf'it) show a very similar behavior. The total energy of both wave packets is 
24.76 GeV for our numerical choice a = 0.1 fm. 
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Fig. 4 displays the corresponding longitudinal field energy wt E \t). Before the wave packets 
start to overlap we find (t) = to a very high precision (10~ 26 ). In the overlap region, 
a hump like structure occurs which grows for increasing values of g. For small values of g, 
Wj^\t) vanishes after the wave packets have passed through each other. After a long time 

( F) 

Wl (t) grows very steeply. Fig. 4 shows clearly that the time between overlap and sudden 
growth shrinks for increasing coupling g. We find, that this time difference scales like 1/g 2 . 

Fig. 4 
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FIG. 4. The longitudinal color electric field energy (t) as a function of time is displayed for dif- 
ferent values of the coupling constant g. The time step width of each time step is At = 03/10 = 0.01 fm. 



Now the question arises whether the energy deposit in longitudinal links can be associated with 
fields propagating into transverse directions. To define energy currents, we denote the Poynting 
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vector in the adjoint representation 



S:=cS x B. 



(25) 



With the transverse and longitudinal components of the vector (25) we define the total trans- 
verse and longitudinal energy currents 



In the following Fig. 5, we show the total transverse energy current ir(t) in the overlap time 
region for different values of the coupling constant. The figure shows that there is no transverse 
energy current before the wave packets start to overlap. In the overlap time region between 
the time steps tsoo and ^1500, we observe a strong increase of the transverse energy current. The 
maximum is reached at time step ti 500 where the overlap is maximal. At decreasing overlap 
in time region between t 150 o and t 2 2oo, ir(t) vanishes. The height of the humps increases for 
increasing values of g and scales like l/g 2 . To become more precise, in the time region of the 
humps as shown in Fig. 5, the transverse current implies an oscillating structure of the time 
period T t = A/4 over which we have averaged to obtain the smooth curves in Fig. 5. A is here 
defined as A = (2n)/k 3 In the upper left panel in Fig. 5, we show an example for k 3 = 7r/(4a 3 ) 
in the case g = 8 in which no averaging has been done. The curve displays a period of T t = 2a^ 
according to A = 8a 3 . No averaging was necessary for times after the overlap time region. 
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(27) 
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Fig. 5 




number of time steps 

FIG. 5. The total transverse energy current ir(t) as a function of time is displayed for different 
values of the coupling constant g. The average longitudinal momentum of the colliding wave packets 
is &3 = 7r/(2ci3). In the upper left window, the full time-dependence of the transverse energy current 
is displayed at the example g = 8 and k% = 7r/(4a3). 

Fig. 6 displays the same transverse energy currents as Fig. 5 but for a larger time interval 
and on a larger scale. Note that the current is plotted in units of [^t(0] = cGeV in Fig. 6 
but in units [«t(0] = cMeV in Fig. 5. After passing through the hump region between time 
step t 80 o and £2000; the energy current txit) disappears. For the smallest coupling g — 4, it 
takes about 3500 time steps or r ~ 35 fm from maximum overlap until to the point where %t 
starts to re-grow. This time, %t grows to much larger values and stays large at later times. The 
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time delay r depends on g : Ak 3: k 3 . We find that r ~ 1/g 2 but haven't studied carefully the 

—3/4 

dependence on Ak 3 and /c 3 . For large A; 3 at least, one can argue analytically that r ~ /c 3 
This will be further discussed in the appendix B. Some few calculations for smaller fc 3 have 
shown that r increases monotonically for decreasing fc 3 . The sudden increase of %t defines the 
"glue burst" . For large coupling it occurs already in the overlap region. 



Fig. 6 
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FIG. 6. Same as in Fig. 5 but for a larger time intervall and on a larger scale for the transverse 

energy current irit)- 

The following Fig. 7 displays the corresponding longitudinal energy currents iL{t). It shows 
that ii decreases at the same time when i T increases. 
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Fig. 7 
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FIG. 7. The longitudinal energy currents iL{t) corresponding to the cases shown in Fig. 6 are 
displayed. iz,(i) decreases at the same time when ixit) increases. 

Subsequently, we present some results from the collision of wave packets with a finite transverse 
extent. We have carried out similar calculations as in the case of plane waves and found a 
similar, but even more pronounced behavior for the time dependence of the transverse and 
longitudinal energy currents. Most of the calculations were carried out on a lattice of the size 
4x200x800 points. The wave packets have been initialized in the same way as in the case for 
plane wave packets. However, the initial wave functions depend now also on the transverse 
coordinates. Since a lattice with the same extension in both transverse directions would lead to 
an exceedingly large number of lattice points, we chose a small number of lattice points in the 
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Xi-direction. Thus, we study the dynamics of a collision in the (x 2 , £3)-plane integrated over 
the Xi-coordinate. Under these restrictions, the scalar function (p(t,x 2 ,x 3 ) defines the initial 
wave packet 

<j>(t, x 2 , x 3 ) := 0o exp (~Ak 2 3 (t + x 3 ) 2 ) exp (-(Ak 2 x 2 ) a ) cos (k 3 (t + x 3 )), (28) 
O := y/(2Ak 3 )/(y/Z<xk 3 ). (29) 

The additional factor exp (— (Ak 2 x 2 ) a ) describes the shape of the wave packet in the transverse 
x 2 -direction. The parameter Ak 2 determines the transverse width and the parameter a controls 
the transverse shape. For large values of a, we obtain a sharp surface and the exponential tails 
are suppressed. Since the transverse extension of the lattice is limited, we choose the extension 
of the initial wave packet in the ^-direction small enough to leave space for the transverse 
dynamics after the collision. Good values are Ak 2 = l/(40a 2 ) and a — 8. The wave packet is 
thus 8 fm wide and extends over 80 lattice points into the main transverse direction. Further 
outside the amplitude is practically zero since a is large. For Ak 3 = 7r/(100a 3 ), the FWHM in 
longitudinal direction is 10 fm. 

With the function 0(t, x 2 ,x 3 ) we determine the initial link variables and the color electric field 
amplitudes and map these data on the lattice at the initial positions of the wave packets. For a 
collision, the wave packets have opposite average momenta ±k 3 . The evolution of the collision 
is carried out in the same way as for plane wave packets. Fig. 8, shows snapshots of the 
transverse energy density of the color electric fields which we define as 

r 2 

w ( T E \t, x 2 , x 3 ) = d Xl J2 Tr(Si(t, x)Si(t, (30) 
J i=i 

The simulation shown in Fig. 8 has been carried out for the non-Abelian case and with a 
coupling g = 6. The other parameter settings are as for plane waves above. The upper left 
picture shows the distribution for the time step t wo where the shape has almost not changed 
as compared to the initial distribution. The wave packets decay into the transverse directions 
while they propagate free over the lattice. The upper right picture displays «4^(£i5oo, x 2, x 3 ) at 
maximum overlap. In the middle left Wj?\t 3500 , x 2 , x 3 ) is displayed. The distribution is almost 
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identical with that obtained in the Abelian case at the same time step £3500- At the time step 
t 40 oo, however, the waves packets start to burst. The last two images correspond to the time 
steps t 50 oo and t 600 o- 

FIG. 8. The transverse color electric field energy density Wj?\t,X2,xs) is shown for six selected 
snapshots taken at the time steps tioo, £1500 > ^3500> ^4000) *5000 and ^000 • The corresponding pictures 
are ordered from the upper left to the lower right. 

Note: The Fig. 8 exceeds the size accepted at the lanl-server and was there- 
fore omitted. We recommend to download the postscript file from the web page 
http://www.phy.duke.edu/~poeschl/ under "Rsearch Related Links" and "Preprints". 

IV. ANALYSIS OF THE YANG MILLS EQUATIONS 

In the following, we present the explanation of the glue burst solution by analyzing the Yang 
Mills equations (1). The details of the discussion are presented in the appendix A. With the 
definition of the field tensor T^ u := V^A V — V v A^ , we rewrite Eq. (1) for the gauge fields in 
the adjoint denotation 

[O„(^/-^i")]_ = (31) 

With the definition of the derivative in Eq. (2) follows 

d^A v = d^A" + igd^A^A^ 
+ ig[A„d»A v }_ 
-ig[A„d»A»}_ 

+ g 2 [A„,[A»,A»}_}_ (32) 

Since the time delayed burstdike behavior of the solution occurs also in the case of colliding 
plane wave packets, we first consider this case. The phenomenon occurs in two steps. It is based 



19 



on a delicate interplay of color degrees of freedom between the color sub-spaces span(Ti, T 2 ) and 
span(T 3 ) and between transverse and longitudinal degrees of freedom in the field amplitudes. 
The first step occurs in the overlap region of the wave packets in both, time and space. The 
Yang Mills equations provide a mechanism which transfers field energy from transverse into 
longitudinal field degrees of freedom. It is explained by the equation for the longitudinal gauge 
field components which follows from Eq. (32) for v — 3. For plane wave packets, only the 
fourth term and the g 2 -term do not vanish on the r.h.s. of Eq. (32) in the case v — 3 and 
therefore the equation 

+ g 2 [A»,[A»,A»}_}_ (33) 

describes the dynamics of the longitudinal fields A 3 during the overlap time. Assuming that 
\A*A%\ < \A^A^\ and |AJA£| < \A*A$\ during the overlap time, we find the following leading 
contributions on the r.h.s. of Eq. (33). 

-ig[A», #U"]_ ~ g(Al&A$ - A*cP A%)T 3 (34) 
g 2 [A {lh [A {1 \A^}_}_ = g^A^A^T 3 - g^A^A^T^ 

+ g'A^Af T 3 - g 2 A 3 (1) A?A^ T 2 (35) 

The above assumption is certainly true for the numerical results presented in the previous 
section which have demonstrated that in the overlap time the longitudinal energy densities are 
smaller by at least one order of magnitude in comparison to the transverse energy densities. 
The indices of the space coordinates in Eq. (35) are now denoted in parentheses in order to 
distinguish from color indices. In leading approximation, the dynamics of the longitudinal fields 
in the overlap time is thus described by the following equation 

+5 2 4)4 1) 4 3) r 3 - ^Vf^r 1 

+ <? 2 4 ) 4 1) 4 3) T 3 - g'A^Af T 2 . (36) 
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The order g term on the r.h.s. of Eq. (36) acts as a source term for the third color component 
of the longitudinal gauge field A 3 (t,x). As soon as the wave packets start to penetrate into 
each other, the space integral over this source term 

J s (f) := / d 3 xg(Al(t,x)d 3 A%(t,x) - A%t, x)d 3 A»(t, x)) (37) 

starts to grow. As Fig. 9 shows in a simplified manner, the derivative d 3 A± is negative in the 
overlap region at the beginning and d 3 A% is positive. Consequently, I 3 (t) > at times before 
the wave packets have reached maximum overlap. After the maximal overlap the derivatives 
change the sign and I 3 (t) < 0. 



Fig. 9 
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FIG. 9. Two colliding Gaussian wave packets of gauge fields A\ and A2 for which 
(f)(t,xs) = 4>oex.p(—l/2Akl(t + X3) 2 ) at the beginning of the overlap. The products A^d^A 1 ^ and 
A^d^Ai (plotted as dashed lines) of the field amplitudes (solid lines) contribute to a source current 
in the longitudinal Yang-Mills equation. 

The Fig. 9 displays in a qualitative manner the situation for the most simple case when A; 3 = 
and the wave packets show no oscillation in position space. The situation shown in Fig. 9 
in principle occurs for each pair of overlapping humps in colliding oscillating wave packets 
for which fc 3 > 0. The corresponding time dependent behavior of the source current 13(f) at 
k 3 = is depicted in Fig. 10. When the penetrating wave packets increase their overlap, the 
longitudinal contribution Iz(t) to the source increases first. Since the derivative on the opposite 
side of both wave packets has opposite sign, the source term changes the sign as soon as the 
two humps have passed the maximum overlap. For k 3 > 0, the wave packets depicted in Fig. 
9 would contain oscillations, i.e. the number of oscillations of the curves in Fig. 10 and in the 
source term would increase accordingly. In this case the envelope of the curves would be of 
gaussian shape rather than the curves itself. 
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Fig. 10 




time [fm] 

FIG. 10. The longitudinal contribution I^(t) to the source is plotted as a function of time (doted 
line). The time dependence of the source term leads to a time dependence of the gauge field ^3 as 
indicated by the dashed curve. The corresponding color electric field is indicated by the dot-dashed 
curve. The contribution from the g 2 -term AAs(t) (solid line) adds to the gauge field resulting in the 
total field As(t) (fat solid line). 

As the source term g(Aj l d 3 A 1 ^ — A 2 ^d 3 Ai) grows in time, the longitudinal gauge field A 3 (t, x) 
grows. At the time of maximal overlap the source current changes the sign and cancels the 
gauge field A 3 (t, x) as shown qualitatively by the dashed curve in Fig. 10. 
During the whole time interval in which the wave packets pass through each other the finite 
longitudinal gauge field A 3 (t,x) enters into the last term on the r.h.s. of Eq. (32) leading to 
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an additional but small contribution to the source term as long as g is not too large. For large 
g 2 , the g 2 -term takes over during the overlap time as soon as the order g term has generated a 
finite amplitude A 3 . 

Our calculations show that for not to large values of g (g < 10) and for small overlap times the 
contribution from the g 2 -term is not dominant compared to the contribution resulting from the 
first source term on the r.h.s of Eq. (36). Furthermore, the g 2 -term can not act as a source 
as long as A 3 = 0. The finite longitudinal amplitude has first to be generated by the fourth 
term on the r.h.s. of Eq. (32) or the first term on the r.h.s of Eq. (36) respectively which 
acts as a initiator for the g 2 -term. The fourth term also again switches off its contribution 
to the longitudinal gauge field when the wave packets have passed the maximum overlap. 
The fourth term goes to zero at vanishing overlap and doesn't contribute any further in the 
receding wave packets. However, the g 2 term causes a finite but small contribution AA 3 (t) to 
the longitudinal gauge field during the overlap time. The qualitative time dependent behavior 
of this contribution is shown by the solid line in Fig. 10. AA 3 (t) which is not canceled after 
the overlap is proportional to the surface under the dashed curve (^4. 3 (t)) in Fig. 10. The fat 
solid curve depicts the time dependent behavior of the total longitudinal field v4. 3 (t) while the 
dashed curve is obtained without g 2 -term. 

As a consequence of this mechanism, a small longitudinal color field component A 3 (t,x) is left 
over in the receding wave packets. This contribution initiates the second step which forms 
the burst itself. Without restriction of the general case, we assume that the first wave packet 
is initially polarized in span(Ti) and the second wave packet is polarized in span(T 2 ). As 
soon as the wave packets obtain a small contribution in the gauge field component A^\ the 
spacial component A^ starts to change the polarization in color space. In wave packet (1) 
the components A^ and A^ grow whereas A^ decreases accordingly. In wave packet (2) the 
components A^ and A^ grow whereas A^ decreases accordingly. This behavior results in 
a rapid change of the T 3 -terms on the r.h.s. of Eq. (36). It is explained by the Yang-Mills 
equation for the first component which follows from Eq. (32) for v — 1. Analyzing the r.h.s. 
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terms, we find 



d^A 1 = -gA^d^AfT 2 + gA^d^AfT 1 

-gAfd {?) A^ + gAfd { , ) A^ + .... 
+ g 2 4 ) A ( i ) A ( l ) T 1 + g 2 AfAfA^T 2 + .... 



(38) 



where we denote only the most important contributions on the r.h.s. . The first four terms on 
the r.h.s. come from the second term on the r.h.s. of Eq. (32). Similar contributions which 
come from the third term are omitted for brevity since they do not bring a new aspect into the 
discussion. The fourth term does not contribute for plane waves. The last two terms on the 
r.h.s. of Eq. (38) come from the g 2 -term. For details, we refer to the appendix A. 
In the following we discuss the mechanism described by Eq. (38). The first, third and sixth 
term on the r.h.s. lead to excitations of modes in the color direction T 2 in the wave packet (1). 
The second, fourth and fifth term act in a similar manner in the wave packet (2) in the color 
direction T 1 , i.e. the amplitude A\ is decreased. The same happens in wave packet (2) but for 

(3) 

exchanged indices. In all six terms appears the longitudinal amplitude A 3 which causes the 
change of the color polarization of the receding wave packets. 

The most important role in this context play the g 2 -terms since the longitudinal gauge field 
enters quadratically. Here, the order g terms on the r.h.s. of Eq (38) act as a trigger for the 
g 2 -term. As for example in the case of wave packet (1), the amplitue A^ is zero initially but 
grows through the corresponding source terms on the r.h.s. of Eq (38). As long as A^ and A^ 
are small, the contribution of the g 2 -term is suppressed by the high powers. If these amplitudes 
exceed a certain strength, the g 2 -term takes over to determine the dynamics of A^ resulting 
in a fast growth. In the case of wave packet (1), these growing amplitudes A^ and A^ enter 
essentially into the third term on the r.h.s. of Eq. (36) inducing the growth of A^\ 
In order to understand the time delay of the burst-like behavior, we study the leading orders 
in the time dependence of the gauge fields shortly after the collision 



4 1) ^4 1) t + ^ 1) t 2 + ... 



(39) 
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Af ^aft + bft 2 + .... (40) 

Both amplitueds are zero before the overlapp of the wave packets. We insert these amplitudes 
into the Yang-Mills equations and consider in leading orders the effects of the g 2 -term on the 



4 3) . 

d^Af * -M^W + - 

= -g 2 a?afa^ + ... (41) 

We integrate both sides of Eq. (41) in time and find a leading fourth order time dependence of 
the longitudinal color electric field = — d t A^ according to 

Ef(t)^W fl | 1 V + -. (42) 

Consequently, the longitudinal field energy increases like Wjf* (t) ~ t 8 in the leading order of the 
time dependence. The high power explains why the longitudinal polarization of the receding 
wave packets is small at short times after the collision but increases rapidly at large times. 
This time-dependent behavior of w[ E \t) together with the resulting transverse expansion of 
the energy densitiy distribution characterises the "glue burst". 

Our numerical results above show that the time-delay of the burst scales essentially like l/g 2 . 
This scaling has the following analytic explanation. Each factor in Eq. (42) is proportional 
to g 2 and thus can be separated as 

<#> = aifV (43) 

where o c W is independent on g. This g 2 -dependence results from the contribution of the g 2 -term 
during the overlap time as discussed above. The equation (42) rereads now 

E?\t) ^ \a?~aM\9 2 ty + (44) 

which explains the numerical observation, i.e. W^ E \t) ~ (g 2 t) 8 . As Fig. 4 shows, the burst 
leads to a peak-like shape in W^ E \t) at the burst time. The steep rise of the longitudinal 
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energy is immediately followed by a strong decrease. This behavior has to be explained by the 
next order in the expansions (39) and (40). The oscillating nature of the solutions suggests 
that the coefficients 6£ have opposite sign as compared to the coefficients a c . This results in a 
next to leading term on the r.h.s. of Eq. (44) which has the form 

+ + + (45) 

Due to the higher power in t, this term takes over shortly after the first term has increased the 
amplitude of E$ (t). The opposite sign however turns the total amplitude back. The following 
drop of the amplitude is stopped by an opposite effekt from the next higher order in Eq. (45) 
and so on. 

In Fig. 4, we observe that w[ E \t) stays finite and seems to oscillate irregularly around a finite 
average value after the burst which we denote preliminary as W L (oo). This average value 
is independent on g which is explained by Eq. (44). Eq. (44) shows that in lowest order in 
time, g 2 can be completely scaled out with t. As discussed above, the lowest order determines 
the rise of the first hump in W^ E \t) at the beginning of the burst. The rise of the first hump 
determines W L (oo). 

The arguments made in the discussion hold for wave packets of large finite transverse extension 
where we can neglect surface effects. For small transverse extensions, the discussion is in 
principle similar but more involved due to contributions from the surface. 

V. SUMMARY AND OUTLOOK 

We have studied time-dependent solutions of the classical Yang-Mills equations which describe 
the collision of initially polarized wave packets in color space and position space. We have sim- 
ulated the collisions on a three dimensional gauge lattice numerically applying the Hamiltonian 
approach of Kogut and Susskind to describe the dynamics of the color fields in SU(2) gauge 
symmetry. As a function of time, we have calculated the transverse and longitudinal energy 
densities w^\t,x 3 ) and w^(t,x 3 ) of the color electric fields. For initially transverse polarized 
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colliding plane wave packtes and for colliding finite wave packets as well, the longitudinal en- 
ergy densities show a strongly time dependent increase in the overlap region around the center 
of collision but vanish when the wave packtes recede. A similar time- dependent behavior was 
found for the transverse total energy current. A certain time r after the collision both the 
longitudinal color electric field energy and the transverse energy current increase rapidly while 
the distribution w^(t,xs) starts simultaneously to decay. Visualizations in three dimensions 
show that the wave packets suddenly decay fast in a decoherent manner when the time r is 
reached. Both, the maximum transverse energy current at total overlap and r scale like 1/g 2 
but for different reasons. This and also the burst can be explained by analyzing the Yang-Mills 
equations. 

The question arises, whether this pure classical phenomenon could play a role in high energy 
nucleus-nucleus collisions. In the present calculations each wave packet was carrying an energy 
of ~ lOGeV. The size of the finite wave packets was 10 fm in longitudinal and 8fm in the 
transverse directions. The energy of ~ 10 GeV is close to the upper limit that can be described 
on a lattice with the constant ai = 0.1 fm for the above size of the wave packet. It is well 
known experimentally that about half of the energy in a nucleus is carried by glue fields. A 
208 P6-nucleus at 100 GeV ~/u carries thus about ~ lOTeV in glue fields. It would be interesting, 
to study the pure glue field dynamics classically in colliding finite wave packets each carrying 
an energy of 10 TeV. This however requires extremely small lattice constants and very large 
numbers of lattice points. The size of the wave packets can be adjusted to the size of colliding 
Pb nuclei, i.e. it should have an extension of about 11 fm in the transverse directions. Such a 
description of course is still very rough and would require much improvement in the future. 
It would also be iteresting to perform a Fourier analysis of the collisions in three dimensions for 
each color separately. There is hope that Dirac-Fermion fields can be included in the future. 
The authors thank S.G. Matinyan for discussions. This work was supported by the U.S. De- 
partment of Energy under Grant No. DE-FG02-96ER40495. 
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APPENDIX A: ANALYSIS OF THE SOURCE TERMS 



We discuss the source terms on the r.h.s. of the the Yang Mills equations which read 

d^A" = d^A" + igd»[A»,A l/ }_ 
+ ig[A ll ,& t A v \_ 
-ig[A^d»A»}_ 

+ g 2 [A„,[A»,A»}_}_. (Al) 

Since the time delayed burst-like behavior of the solution occurs also in the case of colliding 
plane waves, we discuss Eq. (Al) first for this case. The numerical results presented in the 
section 3 have shown that longitudinal energy densities correspond to transverse energy cur- 
rents. The longitudinal energy densities therefore exhibit the basic feature of the glue burst 
solution. The corresponding color electric field components are determined by the negative 
time derivative of the longitudinal components of the gauge fields. We therefore begin with 
the discussion of the time evolution of these field components and we set v = 3 in Eq. (Al). 
Before the collision both wave packets are polarized in the Xi-direction in Euclidean space. 
Accordingly, the longitudinal components are A 3 = before overlap. In the following we focus 
on the time region where both wave packets start to overlap. 
The first term on the r.h.s. of Eq. (Al) contains a sum of four terms 

d^A 1 = d d 3 A° + d 1 d 3 A 1 + d 2 d 3 A 2 + d 3 d 3 A 3 . (A2) 

The first term on the r.h.s. of Eq. (A2) is zero because we use the temporal gauge in which 
A = 0. The second term and the third term are zero because d\A l = and d 2 A 2 = for plane 
waves that propagate into the ^-direction. For v — 3, explicitely considering the contribution 
from Eq. (A2) on the r.h.s., Eq. (Al) now reads 

d d°A 3 + d^ 1 A 3 + d 2 d 2 A 3 + d 3 d 3 A 3 = d 3 d 3 A 3 + ... (A3) 
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This shows that the remaining term of the expression (A2) is canceled by the last term on the 
l.h.s. of Eq. (A3). For plane waves, the second and third term on the l.h.s. are zero. The 
contributions from the other source terms in Eq. (Al) are indicated by the dots. 
The second term on the r.h.s of Eq. (Al) is zero initially because of two reasons. First, A 3 = 
at the very start of the overlapping. The second reason is explained in the following analysis. 
The term in d^A^, A 3 ]_ corresponding to /i — vanishes due to temporal gauge, the /i = 3 
term vanishes because of the commutator, and the /i — 1,2 terms vanish for plane wave packets. 
Consequently, we find 

d,[A^A 3 ]_ = 0. (A4) 
By similar arguments and because of [A3, <9 3 ^4 3 ]_ = 0, we obtain 

iglAp, d»A 3 }_ = ig[Ai, d l A 3 }_ + ig[A 2 , d 2 A 3 }„ + ig[A 3 , d 3 A 3 ]_ = (A5) 
for the third term on the r.h.s. of Eq. (Al). 

The fourth term on the r.h.s. of Eq. (Al) plays an important role in the overlap region. 
According to our calculation for plane waves, we assume that the colliding wave packets are 
initially polarized in color space in the directions Ti and T 2 respectively. When the wave packets 
start to overlap, a superposition of two colors occurs and we obtain 

-ig[A»,d 3 A»\_ = -igUlT, + A\T 2 + ...),9 3 (^T 1 + A»T 2 + ...)]_ 
~ —ig (A^A^ [Ti ,T 2 ]_ + Aj0U?[T 2 , T 1 ],) 
= giAl&AS - Ald 3 A»)T 3 (A6) 

The color indices 1 and 2 refer here at the same time to the contributions of wave packet 1 and 2. 
This assumption is allowed without restriction of the general case in which the wave packets are 
polarized in any direction in the color sub-space span(Ti, T 2 ). On the r.h.s. of Eq. A6, we have 
neglected the terms —ig[A 3 T 3 , d 3 (AiT 1 + A 2 [T 2 )}_ and -ig^A^Tx + AjT 2 ), d 3 A 3 l T 3 }_ since 
A3 = at the begin of the overlapping. Here, for briefness, we assume that (A 3 ^! <^ |A-L4-i | 
and |A 3 Ag I \A 2 A2 1 during the overlap time. Consequently, we neglect excitations of modes in 
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the color sub-space span(Ti,T 2 ) in A 3 (t,x) during the overlap time. These components can be 
included in a straight forward manner in the subsequent discussion but they do not change the 
results qualitatively. What we show does not depend on the above assumption. Nevertheless, 
our numerical results show that the assumption is true for many cases. The remaining term on 
the r.h.s. in Eq. (A6) acts as a source term for the third color component of the longitudinal 
gauge field A 3 (t,x). Its effect is discussed in section 4. 

In the following we analyze the g 2 -term assuming that the remaining longitudinal fields are 
essentially polarized in the T 3 -direction in color space. For briefness, we omit the contributions 
from the first and second color direction. This is sufficient since the T3-component becomes 
always finite in the overlap region as shown above. The g 2 -term breaks up into 

g 2 [A„ [A^A 3 ]_]_ = g 2 [A 1} [A\A 3 ]_]_ + g 2 [A 2 , [A 2 ,A 3 ]_]_ (A7) 

In the following we discuss the contribution of the first term on the r.h.s. of Eq. (7). The 
second term does not contribute if the wave packets are polarized into the Xi-direction or it 
contributes in an analogous manner. In the following, we put indices in Minkowski-space into 
parentheses in order to distinguish from color indices. With the above assumtion, we obtain 

[A\A 3 ]_ ~ [4 1) T 1 + 4 1) T 2 + 4 1) T 3 ,4 3) T 3 ]_ 

= -lA^AfT 2 + lA^AfT 1 (A8) 

for the inner commutator. For Eq. (A7), we thus obtain 

g 2 [A {1) , [^U (3) LL = -g 2 [A {1) ,iA?A^T 2 - iA?A^}_ 

= g 2 A\ l) A i l ) AfT 3 -g 2 A 3 l) A i l ) AfT' 

+ g 2 A 2 {l) A^Af T 3 - g'A^Af T 2 (A9) 

As soon as the wave packets obtain a small contribution in the gauge field component A y 3 , 
the component A^ starts to change the polarization in color space. In wave packet (1) the 
components A^ and A^ grow whereas A^ decreases accordingly. In wave packet (2) the 
components A^ and A^p grow whereas A^ decreases accordingly. 
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This behavior is explained by the Yang-Mills equation for the first component 

d^A 1 = d^A" + igd,[A^A l }_ 
+ tg[A„d^A 1 }_ 
-ig[A„d l A»}_ 

+ g 2 [A»,[A»,A 1 }}_. (MO) 

The first term on the r.h.s. of Eq. (A10) does not contribute for plane waves since d 1 ^ 1 = 0. 
The second term yields the contributions 

igd^A^A 1 ]. = tg[d 3 A^T 3 + ...^T 1 + A^T 2 + ...]_ 
igiA^T 3 + .., dsA^T 1 + d 3 A^T 2 + ...]_ 
= -gA^d^T 2 + gA?d {3) A^ 

-gAfdwdpl* + gAfd^A^T' + .... (All) 

of which the first and third lead to excitations of modes in the color direction T 2 in the wave 
packet (1). The second and fourth term leads to excitations of modes in the color direction T 1 . 
We now discuss the g 2 -term. Here, [^4 1 ,^4 1 ]_ = in the inner commutator and A 2 = for 
polarized wave packets. For the essential contributions of the inner commutator we find 

[A 3 ,A 1 }_ = [d 3 A 3 3) T 3 + ..,4 1) T 1 + 4 1) T 2 + ...]_ 

= iAfA^T 2 - iAfA^T 1 + .... (A12) 

Inserting the r.h.s. of Eq. (A12) into the g 2 -term, we obtain 

AAM' l ,A 1 ]_]_ = g 2 [A i ,[A s ,A 1 ]_]_ 

= g 2 A ( i ) A ( i ) A i l ) T 1 + g 2 AfA ( i ) A ( i ) T 2 + ... . (A13) 



APPENDIX B: TIME DELAY AT LARGE MOMENTA 



In this appendix, we briefly discuss the ^-scaling in the overlap region and in the burst region. 
As discussed in section 4, for not too large g, the fourth term in the equation 
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d^A" = d^A" + igd,\A^A u }_ 
+ ig[A^d»A v }_ 
-ig[A„d»A»}_ 

+ g 2 [A„,[A»,A»}_}_ (Bl) 

acts as the dominant source term for the longitudinal field components in the overlap region 
of the colliding wave packets. We treat Eq. (Bl) first for v = 3. Assuming that during the 
overlap time the wave packets deviate not much from their initial form 

(f)(t,x 3 ) = <f) exp(~Ak 2 3 (t + x 3 ) 2 )cos(k 3 (t + x 3 )), (B2) 

we can make the approximation d 3 A^ ~ k 3 A^ at large k 3 . Consequently, the amplitude A 3 
is proportional to k 3 during the overlap time. This leads to a 1/ /cg-scaling of the hight of the 
hump in (t, z) in the overlap region. Further, A 3 ~ k 3 enters into the g 2 -term on the r.h.s. 
of Eq. (Bl). This leads to a linear /c 3 -dependence of the coefficients which appear in Eq. 
(44) in the section 4. Separating k 3 from the coefficients leads to 

Ei 3) (t) * iitfWWV + - (B3) 
with new coefficients a% = a%/k 3 . We conclude that the time delay of the burst scales like 

v 3/4 . 
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